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We compute single and double inclusive gluon distributions in classical Yang-Mills simulations of 
proton-lead collisions and extract the associated transverse momentum dependent Fourier harmonics 
V2{pt) and vz(pt)- Gluons have a large V2 in the initial state, while odd harmonics such as V3 
vanish identically at the initial time r = O'*'. By the time r < 0 . 4 fm/c final state effects in 
the classical Yang-Mills evolution generate a non-zero V3 and only mildly modify the gluon V2- 
Unlike hydrodynamic flow, these momentum space anisotropies are uncorrelated with the global 
spatial anisotropy of the collision. A principal ingredient for the generation of V2 and V3 in this 
framework is the event-by-event breaking of rotational invariance in domains the size of the inverse 
of the saturation scale Qa. In contrast to our findings in p-|-Pb collisions Yang-Mills simulations 
of lead-lead collisions generate much smaller values of V2,3{pt) and additional collective flow effects 
are needed to explain experimental data. This is because the locally generated anisotropy due to 
the breaking of rotational invariance is depleted with the increase in the number of uncorrelated 
domains. 


I. INTRODUCTION 


A striking result from high-multiplicity proton-proton 
(p-bp) and proton-lead (p-bPb) collisions at the LHC 
was the discovery of a “double-ridge” structure in two- 
particle correlations that is long range in their rapidity 
difference Arj and includes a dominant cos(2A^) modula¬ 
tion in their relative azimuthal angle m- Key aspects 
of the structure of the observed correlations in proton- 
nucleus collisions bear a striking similarity to those ob¬ 
served in heavy-ion collisions [sm and may point to 
a form of collective behavior where many particles are 
correlated with each other. 

In heavy-ion collisions the azimuthal structure of 
multi-particle correlations is quantitatively well de¬ 
scribed by viscous fluid dynamic calculations with fluctu¬ 
ating initial state geometries [12]. This naturally leads to 
the assumption that the physics responsible for the ridge 
structure in high multiplicity p+p and p-fPb events may 
also be driven by the final state collective flow of the sys¬ 
tem. Indeed, calculations using hydrodynamics (or mi¬ 
croscopic models of final state interactions) are able to 
reproduce features of the azimuthal structure in p-|-Pb 
collisions [IMII- Specifically, the Fourier coefficients 
V 2 {pt) and v^ipp) in the expansion of the particle distri¬ 
bution 

dl^ 

^— cx 1 -|- 2 n{pT)cos[n{(l)-i^nipr))] , (1) 

^ n—1 


can be described by several of these models, albeit with 
different parameterizations of the initial state and trans¬ 


port coefficients. In Eq. (1), ipnipr) = ^ arctan 


(sin(n0)).j, 

(cos(n0))0 


is the (transverse momentum dependent) event plane an¬ 
gle associated with the n*** harmonic, and {■)^ denotes 
the azimuthal average with respect to the single inclu¬ 
sive particle distribution at a given pp- 

While some detailed features of the data - such as the 
dependence of V 2 {pp) on the particle mass - can be ex¬ 
plained quite naturally by final state collective effects, a 


number of conceptual problems can be identified within 
this theoretical approach. One concerns the applicability 
of viscous hydrodynamics due to large pressure gradi¬ 
ents na El] that are present for a significant fraction 
of the space-time evolution. Another concerns the sen¬ 
sitivity to the initial state in small systems [20]; which 
requires a better theoretical understanding of the early 
time dynamics. Further, the observation of pronounced 
azimuthal anisotropies even at high transverse momenta 
Pt ^ 3 GeV challenges the hydrodynamic paradigm 
which is best applied to a description of low momentum 
excitations. It is therefore important to explore if multi¬ 
particle correlations at different transverse momentum 
scales can be understood in part or whole in alternative 
approaches. 

For instance, computations of intrinsic two particle 
correlations in the Color Glass Condensate (CGC) frame¬ 
work have been shown to produce azimuthal anisotropies 
compatible with ridge data for even harmonics for pp > I 
GeV in p-l-p and p-bPb collision systems without the need 
for additional final state collective effects [2IH3S]. How¬ 
ever, no odd harmonics were generated because rescat¬ 
tering contributions to the intrinsic correlations were 
not included |26j . In addition to these intrinsic two- 
particle correlations, the presence of domains of directed 
chromo-electric fields inside the proton and the nucleus 
breaks rotational invariance on an event-by-event basis 
and thereby generates azimuthal anisotropies [2Z]. For 
quarks scattering off a colored target, these colored do¬ 
mains generate both V 2 and V 3 [251150] . However, gluons 
scattering off this target only generate even harmonics. 
We note that the possibility that azimuthal anisotropies 
could arise from the event-by-event breaking of rotational 
invariance by color fields has also been considered in a re¬ 
lated approach |5T] . 

Thus far the studies of initial state correlations gen¬ 
erated by color-electric domains in [25] |2S] , and recent 
extensions thereof in |32| . have been based on a descrip¬ 
tion of the proton as a dilute projectile of valence quarks 
scattering off the nucleus. However, experimentally sig- 
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nificant anisotropies are only measured in events with 
very high multiplicities, where it is more appropriate to 
describe both the proton and the lead nucleus as dense 
colored objects. Because of the high gluon occupancies 
in both the proton and the lead nucleus, they can be 
approximated as classical gluon fields, and their leading 
order space-time evolution is described by solving classi¬ 
cal Yang-Mills equations. 

Such classical Yang-Mills (CYM) simulations were per¬ 
formed previously for nucleus-nucleus collisions [ 331131 ] 
(including an early study of elliptic flow in [35]) and 
proton-nucleus collisions [31 [37]. More recently, these 
studies were extended to include more realistic initial 
conditions in the IP-Glasma model [31 [31, which pro¬ 
vides a satisfactory description of multiplicities in high 
energy hadron collisions [37j. When combined with a hy¬ 
drodynamic evolution (MUSIC), the IP-Glasma-l-MUSIC 
model has also been successfully applied to the descrip¬ 
tion of azimuthal harmonics in nucleus-nucleus colli¬ 
sions m-: azimuthal anisotropies in p-fPb collisions were 
however largely underestimated |20] . 

While the azimuthal anisotropy in the aforementioned 
studies [Ullil was generated via the hydrodynamic evo¬ 
lution of the system, the IP-Glasma model also includes 
fluctuations of color charges inside the projectile and tar¬ 
get. Such fluctuations break rotational invariance on 
an event-by-event basis and lead to a momentum space 
anisotropy already present in the initial state. In this let¬ 
ter, we will study these azimuthal anisotropies of gluons 
in the initial state and during the early time dynamics of 
proton-nucleus and nucleus-nucleus collisions. While the 
effect in Pb-fPb collisions is small, we find that initial 
state and early time effects are sizable in p-fPb collisions 
and should be taken into account in the theoretical de¬ 
scription of small collision systems. 

This letter is organized as follows. In Sec. II we briefly 
outline the theoretical framework underlying the classi¬ 
cal Yang-Mills simulations. Our discussion follows the 
literature in the context of the IP-Glasma model [33] 
with the significant modification that we will also con¬ 
sider ‘eccentric’ proton configurations following [40]. We 
then discuss the measurement of azimuthal anisotropies 
in this framework in Sec. Ill and present numerical re¬ 
sults for azimuthal Fourier harmonics of gluons in proton- 
nucleus collisions in Sec. IV. We investigate the sensitiv¬ 
ity of our results with respect to variations in the spatial 
color structure of the proton and perform a comparison 
of the effects in proton-nucleus and nucleus-nucleus col¬ 
lisions in Sec. IV. The final section summarizes our con¬ 
clusions and their implications for collective dynamics in 
proton-nucleus and nucleus-nucleus collisions. 


II. THEORETICAL FRAMEWORK 

Within the CGG framework [3TI143| . the dynamics of 
a high-energy collision is - to leading order in as - de¬ 


scribed by solutions of the classical Yang-Mills equations, 

[D^,F^‘'] = r, (2) 

in the presence of an eikonal color current J''. Here, 
is the covariant derivative in the presence of the field 
and is the gluon field strength tensor. For a 
right moving (projectile) proton and left moving (target) 
nucleus one has 

J'' = (5'"+pp(a;“,xj_) -h (5''“ppb(a:+,xj_), (3) 

and each event is characterized by a different color neu¬ 
tral distribution of random color charges Pp/pb(a;^,x) 
inside the proton and nucleus. 

Before the collision, the small-x gluon fields inside the 
target and projectile nucleus are determined by the solu¬ 
tion of (|^ and can be compactly expressec0 as [441146] 

^p/Pb(xj_) = ^Vp/pb(xj_)5Wj/p|^(xj.) , (4) 

in terms of the fundamental Wilson lines Up/pb(xj^) of 
the projectile and target nucleus. By dividing the lon¬ 
gitudinal direction into Ny discrete rapidity intervals, 
these can be computed as m 

^p/Pb(x±) = nexp('-fg ^^f^^^^^ ^ , (5) 

2=1 X _L / 

for a given configuration of color chargesHere, = 
We employ a McLerran-Venugopalan type model 
[44] for the color charge densities, which follow local 
Gaussian distributions with variance 

g^{p‘^{Yi,x)p\Yj,y)) = S'p/pb(b) ^(^^(x-y) , 

( 6 ) 

where b = (x -|- y)/2. The spatial distributions of color 
charge inside the proton and the nucleus are described by 
S'p/pb(b) and the respective models are outlined below. 

A. The proton 

We will consider two different models for the distribu¬ 
tion of color charge 5'p(b) inside the proton to study the 
effect of the spatial sub-structure of the proton on the 
observed correlations. 

The spherical proton model is a variant of the IP- 
Sat model [33] where the color charge distribution inside 
the proton is spherically symmetric in impact parameter 


^ These expression are valid in light-cone gauge. 

^ We have introduced an effective mass m to regulate the non- 
perturbative large distance behavior. If not stated otherwise we 
will use a fixed value of m = 0.4 GeV in the following. 
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space. In this case, the color charge density is propor¬ 
tional to the saturation scale, i.e., 

S'p(b) = c X t{h)) . (7) 

Qs itself depends on the transverse position b via the 
nucleon thickness function 

The Gaussian width is related to the (two dimensional) 
proton radius relevant to strong interactions as Rp = 
\/2 ■ Bq US]. Its value Be = 4 ± 0.4 GeV“^ as well 
as the dependence of on the nucleon thickness t{h) 
and the center of mass energy ^/s have been extracted 
from fits to deep-inelastic scattering (DIS) data [50] We 
set Be — 4GeV“^ and the only free parameter in this 
model is the proportionality factor c in Eq. Q, defined 
as where is the color charge square per 

unit transverse area. We choose to use a fixed value of 
c = 2 and shall comment later on the sensitivity of our 
results under variation of c. 

The constituent quark proton model was previously 
outlined in [JD] , where the distribution of the color charge 
density is concentrated around the (transverse) positions 
xcQ of three constituent quarks 


where the thickness function T(b) of the nucleus 

A 

T(b) = ^t(b-x,) (II) 

i=l 

is the sum of thickness functions t(b — x^) of individual 
nucleons. We note that for the case of a single nucleon 
{A = 1), this reduces to the spherical proton model and 
we employ precisely the same parametrization in both 
cases. 


C. Early-time dynamics and gluon distribution 

Solving the classical Yang-Mills equations outside the 
forward light-cone leads to the initial state immediately 
after the collision (r = 0+). The initial gauge fields in 
Fock-Schwinger gauge A'^ = A~ + x~ A^')It = 0 are 
given in terms of the projectile and target fields as [511 [51] 

r=0+ “ ^ ^ 

= I [4i;„7i*Pb] , drAX^,^ = 0 , (13) 

and correspond to longitudinal chromo-electric and 
chromo magnetic fields [53] 


S'p(b) = c X 


3 q" 

27r 


Ncq 

^exp 

n—1 



^CQ j 


, ( 9 ) 


which fluctuate from event to event according to a Gaus¬ 
sian distribution with expectation value (xqq) = Be- 
The gluon distribution around each constituent quark is 
spherically symmetric with a radius denoted by Rcq- We 
will use Rcq = y/Bel2 and adjust the overall strength 
Q = 3 GeV, which yields similar results for the dipole 
scattering amplitude (relevant to DIS) as the spherical 
proton model. 


B. The nucleus 

Since we expect a smaller sensitivity of our results to 
the impact parameter dependent structure of the nucleus, 
we limit ourselves to a single model for the color charge 
distribution inside the nucleus. We first sample the posi¬ 
tions Xi oi A = 208 (for Pb) individual nucleons accord¬ 
ing to a Wood-Saxon distribution with radius and surface 
parameters appropriate for a Pb nucleus. We follow the 
IP-Sat model and set 

Sp^{h) = cxQl{^s,T{h)) , ( 10 ) 


^ In practice one extracts the saturation scale Qs(y/s^t{h)) from 
the IP-Sat parametrization of the dipole scattering amplitude. 
More details on this procedure and the general features of this 
model can be found in [37]. 


B'^\r=o+ = ^xAy - dyA^ - ig[A^,Ay\ . (14) 

Starting from these field configurations, the early time 
dynamics in each event can be determined by numerically 
solving the classical Yang Mills equations in the forward 
light cone. Our numerical implementation is based on 
standard lattice gauge theory techniques and we refer the 
reader to [331 [M] [33] for more details of this procedure. 

We can then extract the gluon distribution at differ¬ 
ent times of the evolution by measuring equal-time cor¬ 
relation functions of the gauge fields. We impose the 
Coulomb gauge condition at the time of each mea¬ 

surement and follow previous works [541155] to compute 
the single particle spectrum by a projection on to trans¬ 
versely polarized gluon modes. The single particle distri¬ 
bution is then given by 


dN 

d'^'kxdy 


1 




(15) 


where = (I, —1, —I, denotes the Bjorken met¬ 

ric, A = 1,2 labels the two transverse polarizations and 
a = 1, • • • , Nl — 1 is the color index. In Coulomb gauge 
the mode functions take the form 







2|k^| 



<^(|k±|r) , 


(16) 

( 17 ) 
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/oN 

where kx = {kx,ky) and Ha denote the Hankel func¬ 
tions of the second type and order a (see [54] for details). 

We note that the above definition of the gluon distri¬ 
bution is such that dA^/d^kx is exactly conserved for a 
non-interacting system. This property is important, as 
it will enable us to clearly distinguish between the prop¬ 
erties of the initial state at r = O’*" and the effect of final 
state interactions at later times. 


III. AZIMUTHAL ANISOTROPIES 


Here A^trig is the number of trigger particles in the mo¬ 
mentum bin under consideration, 


o/ A A 1 dA"™® 

is the signal and 

I dA™“ 

H(A7;,A</.) = rfA7?dA(^ 


( 20 ) 


( 21 ) 


is the background contribution estimated from mixed 
events. 


We will now discuss the measurement of azimuthal 
anisotropies, in particular the extraction of V 2 (pt) and 
v^(Pt) for gluons, using two different methods. One is 
based on a measurement of the single particle anisotropy 
in each event while the other follows more closely experi¬ 
mental measurements based on two-particle correlations. 


A. Single particle anisotropy 


Within the single particle method we determine the 
Fourier coefficients V 2 {pt) and v^ipr) from the azimuthal 
anisotropy of the single particle spectra. Since the lattice 
simulation yields the single particle spectrum at discrete 
values of the transverse momenta kx and ky^ we first per¬ 
form a bi-linear interpolation and divide the data into 
transverse momentum |kx| and azimuthal angle (j) bins. 
While in general both the Fourier coefhcients Vnipr) and 
event plane angles ipn{PT) in Eq. 0 depend on the trans¬ 
verse momentum px and fluctuate event by event, we will 
for the moment disregard the px dependence of the event 
plane (EP) and instead compute the second and third or¬ 
der event plane angles ' 4 ’ 2 {Px^) and tpaiPx^) glnons in 
each event as an average over a reference momentum re¬ 
gion, which is chosen to be 1 GeV < Px^ < 6 GeV as 
discussed below. We then extract the Fourier coefficients 
V 2 {EP){px) and V 3 {EP){px) with respect to the refer¬ 
ence event plane according to 


Vn{EP){px) 


I ^ cos(n(0 - ^n{Px^)) 


■ ( 18 ) 


Since Vn{EP){px) fluctuates from event to event our fi¬ 
nal result is obtained by performing an average of these 
quantities over all events. 


B. Two particle correlations 

The measurement of two-particle correlations in ex¬ 
periments is based on the per trigger yield defined to be 


1 ^ ^(A77,A0) 

Ntrig dApdAcj) ’ B{Ari,A(j)) 


(19) 


Since there are no acceptance and efficiency corrections 
in the theory calculation, we can directly compute the 
quantity [5ilS7] 


27r 


dNP^ 


-^trig-^assoc d^cj) 


-{ki,k 2 ) = 


^(fci,<(.)^(fc2,<^ +A<^)) ’ 

( 22 ) 


where Aassoc is the number of associated particles in the 
momentum bin considered and the event average over the 
product of single particle distributions con¬ 

tains the event-by-event single particle anisotropy as well 
as non-factorizable (connected) two-particle correlations. 
Note that while in experiments a rapidity gap is intro¬ 
duced to suppress jet-like correlations around Ary = 0, 
jet-like correlations are not present in our calculation at 
this order [23] and Ary gaps are therefore unnecessary. 

The Fourier expansion of Eq. (22), 


27r 




-{ki,k2) = 


AtrigAassoc dAcj) 

1+ '^2VnA{ki,k2)cos{nA(j)) (23) 


yields the coefficients VnA, from which we define the two 
particle gluon V 2 {px) and vg(px) to be 

Vn{2PC){px) = __ ^ (24) 

^JVr,Aip^^^,pf) 

as is done by the experimental collaborations to measure 
Vnipx) for inclusive hadrons. We choose the reference 
momentum range for gluons to be 1 GeV < < 6 GeV. 

The upper limit in px of this range extends to somewhat 
larger px than that of experimental measurements for 
inclusive hadrons This is to account very roughly 

for the fragmentation of higher momentum gluons into 
hadrons in the px bin of interest. We will compute 
VnAiPTjPx^) in Px bins of width 0.25 GeV defined sym¬ 
metrically around the quoted px values. We will com¬ 
ment in Sec. V on the sensitivity of our results to the 
reference momentum range. 
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FIG. 1. (Color online) Gluon V2 {pt) (left) and V3 {pt) (right) at different times r = 0.0 — 0.4 fm/c in p+Pb collisions at 
impact parameter 6 = 0 fm in the constituent qnark proton model. Open symbols correspond to the single particle anisotropy 
measurement while filled symbols show the results obtained from two particle correlations. Error bands correspond to statistical 
errors only. Experimental results by the ATLAS and CMS collaboration for inclusive hadrons are also shown as a guideline 
for comparison. 


IV. NUMERICAL RESULTS FOR p+Pb 
COLLISIONS 


We will now present numerical results for the az¬ 
imuthal correlations of gluons in p-|-Pb collisions. We 
first study the behavior for collisions of a constituent 
quark proton with a lead nucleus and focus on colli¬ 
sions with zero impact parameter b = Ofm without any 
further multiplicity selections applied. Our results for 
V 2 {pt) and v^{pt) of gluons - obtained from an average 
over Vevt = 128 events - are shown in the left and right 
panels of Fig. We also show experimental results for 
inclusive hadron V 2 {pt) and vz{pt) from the ATLAS [6] 
and CMS [3] collaborations to provide an estimate of the 
relative magnitude and momentum dependence of the ob¬ 
served correlations. We emphasize that since the gluon to 
hadron conversion is by no means straightforward, we do 
not expect a quantitative description of the data within 
the present framework. 

We find that already at the initial time r = 0+ the 
V 2 {pt) of gluons is quite large and extends up to relatively 
high transverse momenta, which can be seen in both 
the two particle correlation and single particle anisotropy 
measurements. The fact that both methods give rise to 
very similar results points to the fact that the origin of 
the observed V 2 {pt) is due to a breaking of rotational 
symmetry of the single particle spectrum. In other words, 
we conclude that gluons are produced with a preferred 
azimuthal direction in each event. 


While the emergence of a preferred azimuthal direc¬ 
tion is well understood in the context of a collective 
expansion of the system, we emphasize that the non¬ 
zero gluon V 2 {pt) at the initial time can not be at¬ 
tributed to collective flow effects. As a consequence of 
the initial color fields in Eq. (141 the energy momen¬ 


tum tensor at r = 0+ is strictly diagonal T'^i/(xj_) = 
diag(e(xj_), e(x_L), e(xj_), —e(x_L)) and the Poynting vec¬ 
tor E X B characterizing the energy momentum flow van¬ 
ishes identically. Instead, the observed V 2 {pt) should be 
attributed to anisotropic gluon production from the fluc¬ 
tuating color fields inside the projectile and target. 

Our result for V 2 and from two-particle correla¬ 
tions can be separated into three contributions. First, 
the Glasma graph contribution, which corresponds to the 
connected graphs analyzed in IIIH2S1 [MSS]. Second, a 
contribution from disconnected graphs, which takes into 
account event-by-event breaking of rotational symmetry 
similar to the effect discussed in [22H29ll32lll9|, and fi¬ 
nally, a contribution from final state interactions repre¬ 
sented by the Yang-Mills dynamics of the produced gluon 
fields. 

The contribution from disconnected graphs can be un¬ 
derstood in the following simplified picture: One can 
imagine the gluon production process as the partons in¬ 
side the projectile scattering off a high-energy nucleus. 
Each parton experiences the color electric field inside the 
nucleus whereby it receives a momentum kick in a cer¬ 
tain direction. In each event the color electric field of the 
nucleus is characterized by several domains with differ¬ 
ent orientations in azimuthal direction and color space. 
While partons scattering off different domains receive a 
different kick, partons (with the same color charge) scat¬ 
tering off the same domain will receive a kick in the same 
azimuthal direction. Hence, depending on the number of 
uncorrelated domains probed by the projectile, the initial 
state particle production can be significantly anisotropic. 

The Glasma graph contribution directly probes the 
color structure within each domain and corresponds to 
genuine non-factorizable two-particle correlations. De¬ 
pending on the degree of polarization of a single domain, 
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the two contributions can be of comparable size [59) . 

With regard to the effect of final state interactions 
on the observed gluon V 2 {pt), we find that the classi¬ 
cal Yang-Mills evolution leads to slight decrease of the 
observed V 2 {pt) with time. This decrease may be at¬ 
tributed to gluon rescattering after the collision. The 
overall effect of the classical Yang-Mills evolution on V 2 
is very small and the initial state value provides a very 
good estimate for the correlation at later times. We note 
again that this is conceptually quite different from a hy¬ 
drodynamic mechanism, where V 2 (jpt) is gradually build 
up as a function of time. 

We now turn to a discussion of the gluon v^(pt) shown 
in the right panel of Fig. [l]for which the interpretation is 
drastically different. We find that the initial state V 3 {pt) 
vanishes identically at r = O'*". This is a consequence 
of the vanishing transverse color electric (Ej^) and color 
magnetic fields (Bj^) at the initial time. With only the 
longitudinal color electric {E'^) and color magnetic (S’*) 
fields being non-zero at r = 0+, the gluon spectrum in 
Eq. (15) takes the form 


dN 


d'^'k^dy 


1 




-0* or s 


7r|k_Lp 


-b 


E2{k^)E2{-k, 


7r|kx| 


(25) 


The above spectrum is manifestly symmetric under ky 0 
—ky “ hence odd Fourier harmonics vanish at the initial 
time. As a result, we can attribute the observation of 
a non-zero vs^px) in our framework at later times ex¬ 
clusively to coherent final state effects included via the 
classical Yang-Mills evolution. 

Quantitatively, we find that already at very early times 
T ~ 0.2 fm/c the classical Yang-Mills evolution leads to 
the build up of a sizable gluon v^(pt) extending to rela¬ 
tively large transverse momenta. Beyond r = 0.2 fm/c, 
the U3 of hard gluons remains approximately constant 
while the low momentum U3 continues to show an in¬ 
crease with time. When we follow the classical Yang- 
Mills dynamics to even later times the system becomes 
more and more dilute and approaches a free streaming 
behavior around r = 0.4 fm/c as previously reported in 

m- 

We note the agreement between the two particle cor¬ 
relation and single particle anisotropy measurement of 
v^(Pt) again points to a correlation between many parti¬ 
cles in each event. While it may appear suggestive that 
the build up of energy-momentum flow may cause a non¬ 
zero V 3 {pt) at later times, to our surprise, we did not 
observe a significant correlation between the global ini¬ 
tial state eccentricity £3 and the px integrated momentum 
space anisotropy V 3 on an event-by-event basis. However, 
it is possible that there is a correlation between the final 
V 3 and triangular anisotropies on shorter length scales 
than the ones probed by the global £3 - such geometri¬ 
cal features are however difficult to extract. Generally, a 
simple description of the non-linear dynamics underlying 


the emergence of a non-zero V 3 {px) remains elusive - a 
deeper understanding of this effect is clearly desirable. 


V. SENSITIVITY TO PROTON STRUCTURE 
AND COLLISION GEOMETRY 

We will now study the effect of the collision geometry 
and system size on the azimuthal correlations. Our re¬ 
sults are compactly summarized in Fig. which shows 
a comparison of the gluon V 2 (px) and v^ipx) at time 
T = 0.4 fm/c after the collision. 

We first analyze the effect of a more detailed sub¬ 
structure of the proton projectile on the correlations in 
p-bPb collisions by comparing a spherical proton with 
one composed of three valence quarks. Generally, a finer 
substructure leads to larger V 2 and V 3 at transverse mo¬ 
menta pt ^ 2 GeV - corresponding to wave lengths on 
size scales much smaller than the nucleon size. How¬ 
ever the overall effect of the proton’s geometry for the 
observed azimuthal anisotropy of gluons is far less signif¬ 
icant than one would expect in a mechanism that gen¬ 
erates azimuthal anisotropies via final state collective ef¬ 
fects [5D]. Since the origin of the observed correlations 
is due to the microscopic structure of color fields, one 
instead expects the correlations to be approximately in¬ 
dependent of the global event geometry. Our result in 
Fig.[^ confirms this picture. 

We have also considered variations of the non- 
perturbative mass scale m and the coefficient c in Eq. Q 
by a factor of two. While in both cases we did not ob¬ 
serve a significant effect on the overall magnitude of the 
observed correlations, we found that for smaller (larger) 
values of m the correlations extend over a slightly larger 
(smaller) transverse momentum range. We note further 
that changing the reference momentum range can 
also have a significant effect on the transverse momen¬ 
tum dependence of the signal which is somewhat more 
pronounced for V 2 (px) as compared to v^^px). While 
the gluon spectrum is generally anisotropic up to high 
momenta, the event plane angles at different transverse 
momenta ipnipx) are not strongly correlated across the 
entire range of momenta. Considering a larger (smaller) 
reference momentum range can therefore lead to a slower 
(faster) decrease at high momenta. Similar observations 
have recently been reported in a related study |32j in 
which the proton was treated as a dilute object. 

We conclude our study with a comparison between 
proton-nucleus (p-|-Pb) and nucleus-nucleus (Pb-|-Pb) 
collisions. In the latter case, we have analyzed collisions 
at two different impact parameters corresponding to cen¬ 
tral (6 = 0 fm) and peripheral (6 = 11 fm) collision|^ 


r When classified in terms of centrality percentile these correspond 
approximately to the 0 — 5% and respectively 50 — 60% centrality 
classes. 
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FIG. 2. (Color online) Dependence of gluon V 2 {pt) (left) and vz{pt) (right) at time r = 0.4 fm/c on the collision geometry 
and system size. Shown are results for p+Pb collisions with two different models for the proton structure and Pb+Pb collisions 
at two different impact parameters. 


We find that in both central and peripheral nucleus- 
nucleus collisions the correlations between gluons are 
much smaller compared to p-|-Pb collisions. Qualitatively 
this difference can be understood when considering the 
different number of localized domains of fluctuating color 
fields responsible for the production of gluons in the dif¬ 
ferent collision systems. While in proton-lead collisions, 
particles are produced from a small number of different 
domains inside the nucleus, the larger overlap area in 
lead-lead collisions gives rise to particle production from 
a much larger number of different domains. Since dif¬ 
ferent domains are uncorrelated with each other the az¬ 
imuthal anisotropy of the gluon spectrum decreases with 
the number of domains (see also m)- Consequently, the 
initial state momentum space anisotropy is much smaller 
in Pb-|-Pb collisions as compared to p-|-Pb collisions. 


VI. CONCLUSIONS 

We have presented results for the azimuthal anisotropy 
of the single and double inclusive gluon distributions 
in p-fPb and Pb-|-Pb collisions obtained from classical 
Yang-Mills simulations. Both the proton and the nucleus 
have been treated as dense QCD objects with high gluon 
occupancy. This description is appropriate for the early 
time space-time evolution of high multiplicity p-|-Pb as 
well as heavy ion collisions at high energies. 

Gluons produced in p-fPb collisions show a significant 
V 2 {pt) already at the initial time immediately after the 
collision. Further evolution governed by the Yang-Mills 
equations modifies this V 2 (pt) only slightly. In contrast, 
odd harmonics of gluons are initially exactly zero, but 
significant values of v^{pt) are built up within times 
T ^ 0.4 fm of the classical Yang-Mills evolution. These 
momentum space anisotropies at early times are uncor¬ 
related with the global spatial anisotropy, in contrast to 


anisotropies generated by collective flow. 

Our results indicate that in p-|-Pb collisions there are 
significant contributions to V 2 and vs from the initial pro¬ 
duction (in the case of V 2 alone) and the early time non¬ 
equilibrium dynamics within the first half fermi of evolu¬ 
tion. These effects cannot be neglected and any calcula¬ 
tion based merely on final state effects is thus incomplete. 

A similar analysis of Pb-|-Pb collisions reveals a dif¬ 
ferent picture. Initial and early-time contributions to V 2 
and V 3 are small, indicating that for larger systems final 
state collective effects are indeed the dominant mecha¬ 
nism for generating the observed anisotropies - at least 
at momenta Pt < 2 GeV, where the presence of such ef¬ 
fects is very plausible. The difference between p-|-Pb and 
Pb-l-Pb collisions can be understood as a consequence, 
in this framework, of anisotropies being generated due 
to localized domains of color fields. A large number of 
mutually uncorrelated domains probed in Pb-|-Pb colli¬ 
sions leads to a nearly isotropic gluon spectrum. Hence 
initial state contributions to V 2 and V 3 are small for large 
collision systems. 

While our study provides a first attempt to quantify 
the importance of initial state effects in high-multiplicity 
proton-nucleus collisions, we expect systematic compar¬ 
isons of p-l-A collisions with deuteron-nucleus (d-l-A) [5D1 - 
IHH| and ^He-|-Au [HI] collisions at RHIC to provide fur¬ 
ther insight into the relative significance of initial state 
and final state effects in small systems. Forthcoming 
p-|-A collisions at RHIC will also help clarify the role 
of nucleon fluctuations relative to sub-nucleon scale ef¬ 
fects in small systems besides providing a comparison of 
results for identical systems at vastly different center of 
mass energies. 
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